Spin-glass-like behavior caused by Mn-rich Mn(Ga)As nanoclusters in GaAs 
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We simulate the indirect exchange interaction between Mn-rich Mn(Ga)As nanoclusters in GaAs 
by analytical means. In contrast to the conventional Ruderman-Kittel-Kasuya-Yosida (RKKY) for- 
mula which considers the mediation by the carriers in the medium, we also include the contribution 
from those inside the clusters. Since the carrier concentration is higher in the clusters, this mod- 
ification allows the RKKY oscillation to change sign. Consequently, while the previous approach 
only favors ferromagnetism for this system, an antiferromagnetic coupling is in fact possible. Since 
the Mn-rich Mn(Ga)As nanoclusters are naturally formed and bound to have different sizes, their 
spin orientation is likely to be frustrated due to mixed preferences from different neighbors. We 
argue that this is likely the source of the spin-glass-like behavior which plagues this system. By 
tuning the size and narrowing its distribution, normal ferromagnetism can be restored with an Curie 
temperature higher than previously thought. 

PACS numbers: 68.35.Ct, 68.35.Fx, 68.35. Rh, 46.65.-|-g 
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DMS (Diluted Magnetic Semiconductors) have been 
hailed as a potential spintronic device to integrate the 
computing power of CPU and the storage ability of the 
hard disc [ll, [3, [I]- By intimately coupling them, we hope 
to come up with a device that can compute while record- 
ing, and vice versa. Besides the ion-implanted samples, 
nanoclusters natually formed or embedded in a semicon- 
ducting matrix have attracted increasing interest due to 
their potential of raising the Curie temperature and rich 
magnetic properties 0, [^, Q, e.g., an enhanced magneto- 
optical and magnetotransport response that can be ap- 
plied to the magneto-optical spatial light modulators 
for volumetric recording. When nanoparticles are small 
enough to enter the superparamagnetic regime, they be- 
have like a giant spin without a permanent magnetic mo- 
ment. 

The indirect coupling between these magnetic 
nanoparticles has been studied in the context of GMR 
(Giant Magnetoresistance) for a metallic medium. 
Note that, although the underlying mechanism is ex- 
pected to follow the theory by Ruderman, Kittel, Kasuya, 
and Yosida (RKKY), how the actual coupling between 
nanoparticles differs from the conventional one, namely, 
between point impurities is not clear. Previous workers 
assumed the validity of conventional RKKY formula and 
concentrated on the mathematical difficulty of integrat- 
ing over spins in a sphere wire, slab, or semi- infinite 
plane [lo|. Discrepancies [ill were found between their re- 
sults and the experiments^, [111, [H, [13] . We believe the 
root of this inconsistency lies in their failure to acknowl- 
edge the mediating role of the carriers inside nanoparti- 
cles. 

This neglect is particularly fatal when the medium 
is semiconducting because, while the previous approach 
would assume a ferromagnetic coupling for all pairs of 
spins, the revised formula actually allows the possibil- 
ity of being antiferromagnetic if only the pair distance 
varies by as little as a quarter of 5 '--^ lOA - typical Fermi 
wavelength for metals. We shall focus on the realistic sys- 



tem of Mn-rich Mn(Ga)As nanoclusters in GaAs[ll, [Tl| 
in this work, and argue that this oscillatory feature is 
intimately related to the spin-glass-like behavior often 
observed in this kind of systcm[l^, [l^ besides the al- 
ready been proposed temperature effect p7j. Armed by 
this knowledge, we explore the prospect of using the nan- 
oclusters to achieve a higher Curie temperature. 
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FIG. 1: Band structure inside the cluster is split by the Zee- 
man energy 2Vo due to the internal field, while that in the 
matrix is assumed to remain a single band. 

When the long-range magnetic order is established 
among nanoclusters, we have to consider the Zeeman 
splitting for the charge carriers. Since the internal field 
in Mn(Ga)As is known to convert its band structure into 
half-metallic [l^l, we duly assign two bands for the car- 
riers inside the cluster, as shown schematically in FiglU 
and different Fermi momenta kp„ can be defined for each 
spin (T. In contrast, the splitting in the medium is rela- 
tively weak and can be neg lected[ll. Let us first calcu- 
late the indirect coupling between a spherical magnetic 
nanoparticle and a single impurity spin. The general 
form of the electron eigenfunction is easiest expressed 
in spherical coordinates with the origin set at its center. 
Symmetry allows us to separate the variables into 

^lrn.= Sl^ir)Yirni0,<j)) (1) 

where the radial part consists of the spherical Bessel func- 
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tions, ji and n/. Since the latter diverges at the origin, 
we have to exclude it for r < a: 



'"^ ' \Ai^ji{k'r) + Bi^ni{k'r), ii r > a 
where the three momenta are related by 
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and the coefRcients can be determined by matching the 
boundary condition as: 



Ala = (k'af 
Bi, = {k'af 



-jjji+i{kaa)ni{k'a) - ji(kaa)ni+i{k'a) 
k 

ii{kaa)ji+i{k'a) ~ -^ji+i{k„a)ji{k' a) 



Finally, after properly normalizing the eigenfunction 
in Eq.© by = I / ^J Af^ + B'f^ , we denote it by 

"^k'lmp- Spatial variation of the magnetization can be 
calculated from the eigenfunction obtained above: 



(^^(^)) - f E E 



n'lml'^k'lrnll 



k' <k' Ini 



(3) 

where s is the spin of carriers and the expectation value 
is calculated in normal ordering. One fact that comes in 
handy is that these nanoparticles shall be sparsely spaced 
in the matrix, namely, their seperation R is much greater 
than their size a. This allows us to take the asymptotic 
limit of the special function by Taylor expand and retain 
just the lowest order term in a/R: 

(MiR)) . 2,sfe-3yQlTl^ '-»-(^^^^^^^+ f ^)) (4) 

where (j> = — arctan((5/VF) and Q{k'p) and W{k'p) are 
constants that depend on the parameters (fc^, a, Vo)'. 



Qik'p) 
W{k'p) 



-A*i^A^+ Bi^ -A*i^Aii+ B^^ Bi 



Nil 

AnBii+Ai,B*^ 



Nil 



Nil 



(5) 



Eq. ([4]) resembles the usual RKKY formula and is sim- 
ilar to the previous result [1, |^, [l^, US] which integrated 
over all pairs of spins but did not include the contribu- 
tion of the carriers inside the clusters. While it is not 
surprising to retain the oscillation with respect to 2kpR, 
we find the energy obtained from Eg. (pi) only becomes 
identical to that of the previous work[8|, 19|, l20| when we 
turn off the Zeeman splitting inside the nanoparticles and 
set its dispersion relation to be the same as that in the 
medium. This implies that (1) our action of matching 
the boundary condition is equivalent to their carrying 
out the integration over spins inside the clusters; (2) the 



coefficient, 2\/Q2 + VF^, in Eq.© can be thought of as 
an effective magnetic moment for the cluster; (3) what- 
ever new physics that may come out of our inclusion of 
the inner carriers must lie in the phase term (p. 

By comparing Eq.(|3]) with the susceptibility expres- 
sion, we can write down a Heisenberg-like Hamiltonian 
for clusters by reducing them effectively to point spins: 



H 



E 



2smk'^J§ cos{2k'pR + (pi + (f>j) 



{2k'pRr 



(6) 



where JoSeff = TT^h^\/Q'^ + W'^/mk'p. Detail deriva- 
tions of Eq.® can be found in our previous article [l9j. 

Finally, we are ready to study the realistic system of 
Mn-rich Mn(Ga)As nanoclusters in the GaAs matrix. 
The fact that their interface is well-defined and both 
are zincblende without obvious dislocations lends cre- 
dential to our treating the broundary as being step-like. 
The carrier density and the effective mass in GaAs were 
roughly 10^^ cm^'^ and O.Smg, and Mn density in the 
clusters was about 10%[l5|. When the Mn ions inside the 
nanoclusters are aligned, strength of the Zeeman splitting 
Vo is estimated to be ''NBS x 10%" where N(3 - 3.3 ev 
and S = 5/2 for Mn[l,[i|. 

The interesting phase 4> and the effective moment are 
plotted against the cluster radius in Figl2l The previous 
approach^, [23| which predicted only ferromagnetism for 
this type of DMS irrespective of the cluster size, is incon- 
sistent with the experiment [lH . In contrast, our result 
allows the coupling to oscillate with respect to the cluster 
size because the Fermi wavelength in the clusters is much 
shorter than in the matrix. Resonance-like peaks arise 
in Figl^Jb) and (c) whenever cluster radius satisfied the 
constraint of quantum interference. Considering the fact 
that the size of these naturally formed Mn(Ga)As nan- 
oclusters is bound to have a distribution, the spin orienta- 
tion of each cluster is likely to be frustrated due to mixed 
preferences from different neighbors. We believe that this 
frustration may be responsible for the s pin-g lass-like be- 
havior common to this type of svstcm|15l.ll6l|. In order to 
be applicable as DMS, the size distribution has to be nar- 
rowed to surpress the spin-glass-like phase. We propose 
that these nanoparticles be prepared separately and spa- 
tially arranged periodically on a substrate before being 
blanketed by layers of semiconductor. The coupling en- 
ergy of this ideal array of nanoparticles is calculated and 
plotted as a function of the cluster size in FiglJ] It shows 
that the previous approach [1, [2^ is likely to underesti- 
mate the optimal Curie temperature by more than three 
times. Since the embedded nanoclusters in many aspects 
share the same features as the magnetic dopants in DMS, 
it is expected that the system under investigation is also 
fit for a spin injector. 

In conclusion, we calculated the indirect coupling be- 
tween the sparsely-spaced Mn-rich Mn(Ga)As nanoclus- 
ters in GaAs. The conventional RKKY formula between 
point spins is modified by an extra phase term after we 
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FIG. 2: (a): The phase term that appears in the inter-cluster 
coupling is plotted against the cluster radius. Panels (b) 
and (c) show how the effective moment of each cluster varies 
with its size. These properties are calculated at 20nm from 
the cluster center - the average spacing between neighbor- 
ing clusters Dashed line is the moment obtained by the 
first-order Born approximation [l^ which is equivalent to the 
previous approach H, and the solid line is our result. 
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FIG. 3: Dashed line is the coupling energy obtained by the 
first-order Born approximation [10] and the solid line is our 
result. The maximum of the latter can be enhanced to more 
than three times that by the previous approach]^, |^. 



take into account the differnt dispersion relation for carri- 
ers across the boundary of nanoclusters. The spin-glass- 
like behavior observed experimentally can be explained 
by the mixture of ferromagnetic and antiferromagnetic 
coupling due to the nonuniformity of the cluster size, 
as opposed to the inhomogeneous spatial distribution of 
magnetic moments in metal[2l|. If the size distribution 
can be narrowed, DMS composed of magnetic nanoclus- 
ters have the potential of achieving a higher Curie tem- 
perature. 
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